Reconstruction of Black Hole Metric Perturbations from Weyl Curvature 
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Perturbation theory of rotating black holes is usually described in terms of Weyl scalars ^4 and if>o, 
which each satisfy Teukolsky's complex master wave equation and respectively represent outgoing 
and ingoing radiation. On the other hand metric perturbations of a Kerr hole can be described in 
terms of (Hertz-like) potentials $ in outgoing or ingoing radiation gauges. In this paper we relate 
these potentials to what one actually computes in perturbation theory, i.e 1/34, and ipo- We explicitly 
construct these relations in the nonrotating limit, preparatory to devising a corresponding approach 
for building up the perturbed spacetime of a rotating black hole. We discuss the application of our 
procedure to second order perturbation theory and to the study of radiation reaction effects for a 
particle orbiting a massive black hole. 

PACS numbers: 04.25.Nx, 04.30.Db, 04.70.Bw 



I. INTRODUCTION 

The spherically symmetry of a Schwarzschild black 
hole background allows for a multipole decomposition of 
metric perturbations, even in the time domain. These 
were studied originally by Regge and Wheeler [1] for odd- 
parity perturbations and by Zerilli [2] for the even-parity 
case. Moncrief [3] has given a gauge-invariant formula- 
tion of the problem, in terms of the three-metric pertur- 
bations. The two degrees of freedom of the gravitational 
field are described in terms of two waveforms, <fi satis- 
fying a simple wave equation 
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(1.1) 



Here r* = r+2M ln(r/2M- 1), and V e ± (r) are the Zerilli 
and Regge- Wheeler potentials respectively. 

Given the solution to the wave equation (1.1) one can 
reconstruct explicitly both the even and odd parity met- 
ric perturbations of a Schwarzschild background in the 
Regge- Wheeler gauge [4, 5]. This permits a complete de- 
scription of the perturbative spacetime. But the Regge- 
Wheeler gauge is unfortunately not asymptotically flat, 
and in order to extract radiation information for a sec- 
ond order perturbative expansion one has to perform a 
new gauge transformation [6]. Moreover, the desirable 
properties of the the Regge- Wheeler gauge being unique 
and invertible are lost in the case the background is a 
rotating black hole, i.e. represented by the Kerr metric, 
where no effective multipole decomposition is yet known. 

There is an independent formulation of the perturba- 
tion problem derived from the Newman-Penrose formal- 
ism [7] that is valid for perturbations of rotating black 
holes. [8] This formulation fully exploits the null struc- 
ture of black holes to decouple the curvature perturba- 
tion equations into a single wave equation that, in Boyer- 



Lindquist coordinates (t, r, 9, if), can be written as: 
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where M is the mass of the black hole, a its angu- 
lar momentum per unit mass, £ = r 2 + a 2 cos 2 9, and 
A = r 2 — 2M r + a 2 . The source term T is built up from 
the energy-momentum tensor [8]. Gravitational pertur- 
bations, corresponding to s = ±2, are compactly de- 
scribed in terms of contractions of the Weyl tensor with 
a null tetrad. The components of the tetrad (also given 
in Ref. [8]) are appropriately chosen along the repeated 
principal null directions of the background spacetime [see 
Eq. (1.4) below]. The resulting gauge and (infmitesi- 
mally) tetrad invariant components of the Weyl curva- 
ture are given by 
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where an overbar means complex conjugation and p is 
given in Eq. (1.5) below. Asymptotically, the leading 
behavior of the field vb represents either the outgoing ra- 
diative part of the perturbed Weyl tensor, (s = —2), or 
the ingoing radiative part, (s = +2). 

The components of the Boyer-Lindquist null tetrad for 
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the Kerr background are given by 

(H = 



r 2 + a 2 a 



(1.4a) 
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(r 2 +a 2 ,-A,0,a), (1.4b) 



2(r 2 +a 2 cos 2 0) 

(m Q ) = -p (ia sin 0,0, sin 6). (1.4c) 

V 2(r + ia cos 9) 

With this choice of the tetrad the non- vanishing spin co- 
efficients are 

1 n _cot0 
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(1.5) 



and the only non-vanishing Weyl scalar in the back- 
ground is 



V> 2 = Mp 3 



(1.6) 



Analogously to the Zerilli-Regge- Wheeler waveforms, 
ip4 can be directly used to compute energy and momen- 
tum radiated at infinity, but it remains to relate it to 
metric perturbations. Chandrasekhar [9] studied a way 
to to obtain metric perturbations of the Kerr metric, but 
it was proved by Price and Ipser [10] that this choice is 
not a proper gauge, namely is an incomplete constraint 
on the coordinates. 

Chrzanowski [11] generalized work of Cohen and Kege- 
les [12] on Hertz potentials to the gravitational perturba- 
tions of the Kerr metric. In Rcf. [11] explicit expressions 
are given for homogeneous metric perturbations in the 
frequency domain. Wald [13] subsequently showed that 
the expressions given in Ref. [11] do not lead to real met- 
ric perturbations. Cohen and Kegeles [14] then corrected 
their expressions and gave explicit equations (see Sec. II) 
relating metric perturbations to a gravitational Hertz po- 
tential, that fulfills Eq. (1.2), but that is different from 
ip4 or tp . In those works no explicit method was given 
for determining * in any specific astrophysical problem. 

In Sec. Ill we provide the explicit formulae relating 
a gravitational Hertz potential to ip4 or tpo in the time 
domain, hence to the given initial data defining the as- 
trophysical model one wants to evolve (See Ref. [15] for 
the 3+1 decomposition of the Weyl scalars). These allow 
one to define the outgoing and ingoing radiation gauges 
that are asymptotically flat at future infinity and regular 
on the event horizon respectively. Such gauges have been 
found [16] especially well suited for computing second or- 
der perturbations of a Kerr hole and, once generalized in 
presence of matter, can provide a first step toward com- 
puting radiation reaction (self-force) corrections [4] to the 
trajectory of a particle orbiting a rotating black hole. 



Vacuum second order perturbations have a direct ap- 
plication to the close limit expansion of the final merger 
stage of binary black holes with comparable masses. Per- 
turbation theory in conjunction with limited full numer- 
ical simulations has proved to be an extremely powerful 
tool to compute waveforms from binary black holes from 
near the innermost stable circular orbit [17, 18]. The 
combination of radiation reaction correction plus second 
order perturbations provides a formidable tool for com- 
puting gravitational radiation from binary black holc- 
neutron star systems, and is particularly relevant to the 
computation of template banks for ground based inter- 
ferometers such as LIGO/ VIRGO/CEO about to enter 
on-line, as well as space missions such as LISA, sensi- 
tive to supermassive black hole binary systems. Thus in 
order to incorporate these improvements to our theoret- 
ical predictions it is imperative to know how to build up 
explicit metric perturbations around a Kerr background. 



II. FORMULATION OF THE PROBLEM 

We use the notation 

to describe metric perturbations. 



(2.1) 



A. ingoing and outgoing radiation gauges 

Chrzanowski [11], and Cohen and Kegeles [14] found 
two convenient gauges that allow one to invert the met- 
ric perturbations in terms of potentials ^irg ot ^org 
satisfying the same wave equations as the Weyl scalars 
P~ 4 ip4 or tpo respectively. 

In the ingoing radiation gauge (IRG) we have 

h u = h^W = 0; h ln = h^Fn" = 0, 
him = h^m v = 0; h m = h^rr? = 0, 
hmm = h^m^lff = 0, (2.2) 

and the homogeneous (for vacuum) metric components 
can be written, in the time domain, in terms of solutions 
to the wave equation for p~ 4 ip4 only, as follows: [27] 



hnn G - -{((5 + S+3/3-r)( ( 5 + 4/3 + 3r)} 



(^irg) + c.c. 



(2.3a) 



h^ G = -{(D-p)(D + 3p)}(* IRG ), (2.3b) 



h 



IRG 

(nm) 



= -7j{(S -S + 3/3-7f -t)(D + 3p) + 
(D + p-p)(6 + 4(3 + 3T)}(H> IRG ), (2.3c) 



where c.c. stands for the complex conjugate part of the 
whole object to ensure that the metric be real [13, 14]. 
Note that in this gauge the metric potential has the 
properties of being transverse (h^l 11 = 0) and traceless 
(h£ = 0) at the future horizon and past null infinity. It 
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is thus a suitable gauge to study gravitational radiation 
effects near the event horizon. 

The complementary (adjoint) gauge to the ingoing ra- 
diation gauge is the outgoing radiation gauge (ORG), 
which can be obtained from the IRG upon exchange of 
the tetrad vectors I <-» n, m <-> m and the appropriate 
renormalization. It satisfies: 



= hir, 



= h n 



(2.4) 



The metric potential has now the property of being trans- 
verse [h^n^ = 0) and traceless (h£ = 0) at the past 
horizon and future null infinity. It is thus an example of 
a suitable asymptotically flat gauge in which to directly 
compute radiated energy and momenta at infinity. 

In this gauge, the homogeneous metric components can 
be written in terms of solutions to the wave equation for 
ipo, as 



h 



ORG 



= -p- A {(5 -3a- [3 + 5n)(5 - 4a + n)} 



(2.5a) 



ORG 



l ORG 
n (lm) 



4 7 )} 



(^org) + c.c. 

-p- 4 { (A + 5p - 3 7 + 7)(A + n- 
{^org), (2.5b) 
- ^p- A { (5 - 3a + + 5tt + t) (A + ft - 4 7 ) 

+(A + hp - Ji - 3 7 - j)(S - 4a + tt)| 
(*oag), (2.5c) 



where the directional derivatives are D = P<9^,A = 
71^8^,5 = m^d^, and the rest of the Greek letters repre- 
sent the usual notation for spin coefficients. [7] 



B. 4th order equations for the potential 



III. EXPLICIT SOLUTION FOR 
NONROTATING BLACK HOLES 

A. Master equation 

The key observation here is that for a = the differen- 
tial operator C acting on the potentials in Eqs. (2.7) and 
Eqs. (2.9) contains only angular derivatives. Since for the 
spherically symmetric background we can decompose the 
potentials into spherical harmonics of spin weight s and, 
from [19], we have 



{-) m + s _ s Yt.- m , 
s (* + 2)! 



[ -2Yt m ] - (-) m Z ^ - -2Ye- m , 



3333 [ i% m \ 



l ym+2 + 2 )' v 



(3.1a) 
(3.1b) 

(3.1c) 



we obtain a first order relationship among ip4 or ipo and 
the IRG or ORG potentials decomposed into multipoles 



[ 



± 



4{i- 2)! 
2)! 



I *IRG 

3Md t *f RG , (3.2) 



"4(^-2)1 ORG 
where we have used the notation 

^± = I[V/'™±(-)"V 



3Md t *± RG , (3.3) 



(3.4) 



for all fields decomposed into multipoles. 

Since ^irg and ^org satisfy the master equation 
(1.2) for spin s = +2 respectively we can eliminate from 
this equation all time derivatives by replacing Eqs. (3.2) 
or (3.3) and its time derivatives into the Teukolsky equa- 
tion. This leads to the following equation for the IRG / 
ORG potentials, both represented here by \& 



From the evolution of the Teukolsky equation we can 
obtain ip^ and ipo, the two gauge and tetrad invariant 
objects representing outgoing and ingoing radiation re- 
spectively. To relate the unknown potential <i? to them 
we use their definitions (1.3) to obtain 



ipo = DDDD [* JflG ] , and 



(2.6) 



P" 4 V>4 = ^ [CCCC^irg - 12p- 3 ip 2 d t y IRG ] , (2.7) 



where C = <5+iasm9d t and 9 = — [de + s cot9—i csc8d, 
while for the outgoing radiation gauge we have 



ipo 



A 2 AAAA [A 2 * i?c] , and (2.8) 
i [CLCC^orc + 12 p- 3 ^, 2 d t y org] ■ (2.9) 



In order to obtain an expression for the potentials in 
terms of the known quantities ip± or ip , one has to invert 
a fourth order differential equation where ip4 or ipo act as 
source terms. This will be the central task of our paper. 



where 



A- s d r [A s + 1 d r y ± ]--(n AS ) 2 v ± 

A/fr 

±2 S r(04 S )(— -l)tt± 
-(I- s)(l + s + l)$ ± =F ± 



* = 



^ IRG for s = -2 
^ORG for s = +2 



1 (1 + 2)! 
12M {£-2)V 



and the source term is 

„4 



(3.5) 



(3.6) 



(3.7) 



3MA 



3MA 

(Mr/ A - 1)(V' ± ). 



(3.8) 
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The equation (3.5) is our fundamental equation to 
solve for the potential in terms of the known fields ip4 or 
tpo that appear in the source terms. The key observations 
here are to note that, separately for the "plus" and "mi- 
nus" parts: i) the left hand side of this equation contains 
the Teukolsky operator in the frequency domain for the 
algebraically special frequencies, u — ±iflAS, and h) the 
source terms in (3.8) are the Cauchy data for the Teukol- 
sky operator in terms of ip4 or -00, precisely as they would 
appear in a Laplace transform approach to Eq. (1.2) [See 
Eqs. (A2-A3) of Rcf. [20]]. Notably, we have arrived to 
this equation working in the time domain, without any 
frequency decomposition. 



B. solution 

The single frequency appearing on the right hand side 
of Eq. (1.2) is precisely an algebraic special frequency. 
Algebraically special perturbations of black holes excite 
gravitational waves which are either purely ingoing or 
outgoing. Hence only one of ^4 or ip is non-zero while 
the other vanishes. Chandrasekhar [21] has obtained the 
explicit form of the algebraic special perturbations of the 
Kerr black hole. This is a remarkable fact, because we 
can then use these analytic solutions to the Teukolsky 
equation for algebraic special perturbations to construct 
explicit solutions to (3.5). 

For Schwarzschild, Chandrasekhar [21] gives two alge- 
braic special solutions 



2/1 (r) = 
2/2 (r) = 



Xr \ 2 r 2 A 2 (A + l)r 3 

1 + TT + 77TTT7 + 



M 3M 2 



(A+l)r 



9M 3 



3M 



(3.9) 



(3.10) 



where A = {£ — 1)(£ + 2)/2. Note that y 2 gives rise to 
an algebraically special solution for tp4, (s = —2), at 
frequency lu^ s = —iflAS an d for A 2 tp , i s = +2), at 
frequency oj^g = +iflAS, while yi gives rise to an alge- 
braically special solution for tp4 at frequency oj\ s and for 

A 2 V> at frequency lo^s- 

These two solutions satisfy the "plus" and "minus" 
parts of Eq. (3.5) for ^f RG [s = +2) and vice-versa for 
^org ( s = — 2). A second set of independent solutions 
can be found: 

*iW = W(r) f 5S (3.11) 
Jim Vi\ r ) 

<- r W(r') 



z 2 {r) = y 2 (r) 



2/2 (r') 2 

where the Wronskian of the solutions is 



dr' 



(3.12) 



W(r) : = yi, 2 {r)d r z lt2 (r) - z lt2 (r)d r yi, 2 (r) 

= A(r) =r(r-2M). (3.13) 

Note: Z\ and z 2 are not algebraically special solutions, 
although they are each a homogeneous solution to their 



respective Teukolsky equation at an algebraically special 
frequency. 

Making use of these solutions, the explicitly expression 
for the potential can be written as follows: 



*i RG ( r ) = -vi (0 / 

J2 
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OO / / 
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] Sr Mr') 2 

.r^r^^, (3 , 5) 

J r 



+ 2/2 (r) 



Z2(r')F-(r ') ^ 
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valid for each hypersurface where ip or ^4, entering in 
F given by Eq. (3.5), are evaluated. 

The equations for the ORG are obtained by exchanging 
"plus" and "minus" parts above (cf. Eq. (A3) and (A4)), 
and by adopting the corresponding dependence of F on 
the components of ip [See Eq. (1.3)]. 



IV. APPLICATIONS 

A. Metric perturbations (IRG) 

Using the Regge- Wheeler notation [1] for metric per- 
turbations, conditions (2.2) read 

^(r,i)( evon ' odd ) = -(1 - 2M/r) h[{r, t)( cvc ^ odd \ 



H e (r,t) = H e 2 (r,t) = -H((r,t). 



(4.1) 



We can now explicitly compute the metric perturba- 
tions (2.3) or (2.5) in terms of the computed i/jq or ^4: [28] 

[h(nrn)] 1 



{ftS(r,t)(— > - ^M) (odd) } ^J-^ _rY, 
ry x 12M (1 - 2M/r) r 2 J [ IRGl 
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rytJ r A 2[ J 3M(r-2M) 
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+ ry 2 J 2M A 21 J 3M(r-2M)/ 
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x^2(£ + 2)(£-l) _iy /m , 
where C t = (I + 2)1/ (I - 2)!, and 

[h nn ] e = (l-2M/r) H e (r,t) Y £m 

= ^l*iR G }\l¥^ oY lm , (4.3) 
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r,_l(even) _ / , ,, feg] (l2 M 2 + r 2 C;) ^ 

[/lmmJ ~i /b r (r-2M)M 



-1/6 
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(A) (12 M 2 + r 2 C ; ) 
y x r[r -2M)M 



A 2 



[F+] dr' 



r 2 [p- 4 a^ 4 + ] 



2/3- 



.d_\ n -i„l,+ 
dr 



(r-2M) 2 M ~'"(r-2M)M 

1/1t >~V2"] (-36r + 84Af + r 2 Q) 

— l/lo ^ 

(r - 2 M) M 

+ (r 3 C ; 2 + 12 r (6 I 2 + 6 Z + 7 C, - 12) M 2 

-144 (Z + 2) (/ - 1) M 3 - 36 r 2 C\M) 

„ [*k G ] 1 



(r - 2Mf M 2 r 



2 7 
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(-36r + 84M-r 2 Q) 
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(r - 2 M) M 

+ ^ (r 3 C 2 + 12 r (6 1 2 + 6 1 - 7 Ci - 12) M 2 

-144 (J + 2) (/ - 1) M 3 + 36 r 2 C ; M) 



(r - 2M) 2 M 2 r 



(4-6) 



In writing these we have explicitly lowered the order of 
the derivatives of the potential by making use of the fol- 
lowing identities 

d r ([*+]) = 

(^)^-£f(&)^ (4 - 7) 

9 r ([*"]) = 

(^) [*-]+- r (§) m*'. (^) 

V 2/2 / 2/2 7 2M v A 2 / 

and directly from Eq. (3.2) 



t IRG 



l2M(£-2)\ 



3Mp' 



-V>4- (4-9) 



B. Metric perturbations (ORG) 

Using the Regge- Wheeler notation [1] for metric per- 
turbations conditions (2.4) read 

/^M)(° von < odd ) = (1 - 2M/r) h{(r,t)^ ven ' odd \ 



'+!) 

H*(r,t) = H}(r,t) = H{{r,t). 



(4.10) 



The explicit metric perturbations can be found directly 
from the previous subsection, Eqs. (4.2)-(4.6) upon ex- 
changes of the tetrad contractions I — ► n and m — > m, 
and the consequent change in normalizations, s, and po- 
tentials, as described throughout the paper. 



V. DISCUSSION 

We have explicitly computed the metric perturbations 
of a nonrotating black hole in terms of the Weyl scalars 
ip4 and ipo which can be computed directly by solving 
the Teukolsky equation for any appropriate astrophys- 
ical scenario, given the corresponding initial data. In 
doing so we had to invert Eqs. (2.7) or (2.9). This was 
performed making explicit use of the multipole decom- 
position of the potential, Weyl scalars and metric. The 
extension of this procedure to the rotating background is 
not straightforward, but we have learned some key fea- 
tures: The algebraic special solutions of the Teukolsky 
equation will play a crucial role in finding the solution 
for the Hertz potential in terms of ip^ or ipQ. One can 
see this as follows. In order to invert Eq. (2.7) for the 
potential we first seek out solutions of the homogeneous 
equation. Hence tpi should vanish, but for the solution 
to be nontrivial, V'o, given by Eq. (2.6), must not vanish 
[22]. These two conditions are precisely the conditions 
that define the algebraically special solutions for the po- 
tential satisfying the vacuum Teukolsky equation (1.2). 
These two solutions (in the time domain) should allow 
one to build up the Kernel that inverts the fourth or- 
der equation (2.7). An identical argument applies for the 
ORG potential when working with Eqs. (2.8) and (2.9). 

One main application of this formalism is to go beyond 
first order perturbation theory and compute second order 
perturbations of rotating black holes. In Ref. [16] there 
was developed a formalism for the second order Teukol- 
sky equation that takes the form of the first order wave 
operator acting on the second order piece of the Weyl 
scalar ip^, and a source term build up as a quadratic 
combination of first order perturbations. It is precisely 
to compute this source term that one needs the explicit 
form of the metric perturbations. In Ref. [16] it was found 
that to describe the emitted gravitational radiation the 
ORG gauge is specially well suited. Hence one has to 
solve the the first order Teukolsky equation for ip^ and 
tAq 1 -* , then later to build up the source of the second order 

(2) 

piece of the Weyl scalar ip\ (See Eq. (9) in Ref. [16]). 
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A second important application of the reconstruction 
of metric perturbations around Kerr background is to 
compute the self force of a particle orbiting a massive 
black hole [23, 24] and to compute corrected trajecto- 
ries [4] depending on the perturbed metric and connec- 
tion coefficients along the particle world line. This task 
is left for a future paper. While we know the form of 
the Teukolsky equation in the presence of of perturba- 
tive matter around a Kerr hole (see Eq. (1.2)), we need to 
generalize the equation satisfied by the potential and the 
relationship between this potential and the metric per- 
turbations, i.e. the generalization of Eqs. (2.3) and (2.5). 
In particular we know that not all of conditions (2.2) or 
(2.4) can hold in the presence of matter since they are 
then incompatible with the Einstein equations [25]. 
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APPENDIX A: RELATION BETWEEN IRG/ORG 
POTENTIALS 



We give here a relation expressing the potential ^org 
in terms of the result for ^irg- This was not given in 



[11, 12, 13, 14]. We begin by defining a field x> with spin 
weight —2, through: 

*irg = \ [CCCCx + 12p- 3 ip 2 d tX \ , (Al) 

analogous to Eq. (2.9). The solution for ^org is: 

Vorg = DDDD [x] , (A2) 

also the relation (2.6) between -0 O an d ^irg- A potential 
degeneracy for algebraically special modes was discussed 
briefly in [26]. The explicit solution for the modes \ lm 
can be written, in terms of x ± given by Eq. (3.4), as: 



X + (r) = z 2 (r) \ 

J2 



2M 
oo 



+ 1/2 ( 



X 



l{r) l A(r') 2 
-(r) = -yUr) f 



yi(r')V+(r>) 
A(r') 2 

Z ^ V+ ^ d r> (A3) 



r ziWv-W dr , 



-*<0/ "<•'• (A4) 



A(r') 



cf. Eqs. (3.14) and (3.15) above. Here V ± are given by: 
V ± 



4 4 

^A^fnG)T^(n AS )^f RG ) 



+4^(Mr/A-l)(v^ G ), 
similar to Eq. (3.8), but note the sign differences. 
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